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We consider sound waves in superfluid nucleon-hyperon matter of massive neutron-star cores. We 
calculate and analyze the speeds of sound modes and their damping times due to the shear viscosity 
and non-equilibrium weak processes of particle transformations. For that, we employ the dissipative 
relativistic hydrodynamics of a superfluid nucleon-hyperon mixture, formulated recently (M. E. 
Gusakov and E. M. Kantor, Phys. Rev. D78, 083006 (2008)). We demonstrate that the damping 
times of sound modes calculated using this hydrodynamics and the ordinary (nonsuperfluid) one, 
can differ from each other by several orders of magnitude. 
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I. INTRODUCTION 



In recent years there is a growing interest in studies of neutron-star pulsations. This is related to a number of 
reasons. First of all, the recently discovered high frequency oscillations of electromagnetic radiation during giant 
flares may be associated with the pulsations of neutron stars [l|, [2j. Second, the gravitational- wave detectors, which 
will be able to detect gravitational radiation from isolated pulsating neutron stars, are under construction [H, [j, [5|, [|| . 

For interpretation of the observations, it is important to have a well-developed theory of neutron-star pulsations. The 
construction of such a theory is complicated by the fact that the baryons in neutron-star cores can be in superfluid 
state @, H, H, OH HH, EH- Thus, to study the pulsations, one has to use a hydrodynamics describing mixtures of 
superfluid liquids. There is a substantial body of literature, devoted to pulsations of superfluid neutron stars (see, 
e.g., Refs. [IE [II, [13, Bill All these papers deal with the nucleon npe{p) matter 

composed of neutrons (n), protons (p), and electrons (e) with possible admixture of muons (/i). Most of them study 
pulsations at zero temperature (see, however, Refs. (la. Ilq|). 

In this paper we for the first time investigate dynamic properties of superfluid nucleon-hyperon matter in the cores 
of massive neutron stars, composed, in addition to neutrons, protons, electrons, and muons, of A and S~ hyperons (A 
and £, respectively). For that, we employ the relativistic hydrodynamics [27] . describing a superfluid nucleon-hyperon 
mixture at arbitrary temperature. We study the simplest pulsations in such matter - sound modes, how they travel 
and how they damp. Within this simple example we demonstrate, in particular, that the characteristic damping times 
of pulsations, calculated sclf-consistently in the frame of superfluid hydrodynamics, can differ by several orders of 
magnitude from those calculated using the nonsuperfluid hydrodynamics (in the latter case the effects of superfluidity 
are taken into account only at calculating kinetic coefficients). 

The paper is organized as follows. In Sec. II we give an overview of the main reactions of particle mutual trans- 
formations in the nucleon-hyperon matter. In Sec. Ill we briefly discuss the relativistic dissipative hydrodynamics 
describing superfluid nucleon-hyperon mixture. In Sec. IV we analyze the sound modes in such matter neglecting 
dissipation. In Sec. V we calculate the damping times of sound modes due to the noncquilibrium reactions dH) — (|4j) 
(see Sec. II) and shear viscosity. Section VI presents summary. 

Throughout the paper, unless otherwise stated, we use the system of units in which the Plank constant h, the speed 
of light c and the Boltzmann constant k& equal unity, h = c = k& = 1 • 

II. THE MAIN PROCESSES OF PARTICLE TRANSFORMATIONS IN NUCLEON-HYPERON 

MATTER 



The most effective weak processes in nucleon-hyperon matter are the following nonleptonic reactions [28|, |2j| |3fJ, |31| 

n + n <-> p + E _ , (1) 

n + p <-> p + A, (2) 

n + n «-> n + A, (3) 

n + A <-> A + A. (4) 
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The full thermodynamic equilibrium implies, in particular, the equilibrium with respect to these reactions, 

5(j,i = 2^„ - n P o - fi-zo = 0, (5) 

Sfl 2 = <5^3 = <fy«4 = MnO - MAO = 0. (6) 

In this case the average number of direct reactions in unit volume per unit time is equal to the number of inverse 
reactions. In Eqs. ((5|) and (J6j) fj,io arc the chemical potentials of particle species i = n, p, A, and S, taken at equilibrium; 
S/i m (m = 1,. . ., 4) are the disbalances of the chemical potentials for the reactions In what follows, we mark 

the equilibrium values of thermodynamic quantities with the subscript 0. Accordingly, the thermodynamic quantities 
without the subscript (e.g., /Ltj) refer to perturbed matter. Notice that, the equilibrium conditions for the reactions 
©, ©, and (g]) coincide. 

Eqs. ([5]) and ^ do not hold in the perturbed matter (8fj,\ ^ 0, 8^,2 = <^M3 = £>M4 7^ 0), so that the numbers of direct 
and inverse reactions are not equal. The noncquilibrium reactions (JTJ)— (|4j) act to return the system to the equilibrium 
state. This leads to dissipation of mechanical energy, accumulated in the matter. 

Along with the reactions dH) — d4j) there is a number of weak reactions with leptons. The leptonic reactions (e.g, the 
direct and modified Urea processes with electrons or muons) are much slower in comparison to the reactions (JU)— (J3J) - 
In the interesting range of parameters (temperatures and pulsation frequencies) , they cannot influence substantially 
chemical composition of the perturbed matter and will be neglected in what follows. However, we assume that the 
unperturbed matter satisfies the equilibrium conditions with respect to these reactions, 

(InO = ptpO + MeO, MeO = MpO, (?) 

where fi e o and fi^o are the equilibrium chemical potentials for electrons and muons, respectively. In addition to the 
processes described above, there is a fast nonleptonic reaction due to the strong interaction of baryons, 

n + A^p + Y,'. (8) 

We assume that the perturbed matter is always in equilibrium with respect to this reaction (27l . |30| , 

fytfast = Hn + MA - Mz> - Me = 0. (9) 

It follows from the condition (J9j> , that the chemical potential disbalances for all the four reactions dTJ) — ((4j) coincide, 

8/j, = Sfix = 6(i2 = <fy*3 = ^M4- (10) 

Now the equilibrium conditions ([5]) and ([6]) can be rewritten as 

5fi = 0. (11) 

Below we denote the difference between the average number of direct and inverse reactions (pj,. . .,([4|), occurring in 
the unit volume per unit time, by Ari,. . .,Ar4, respectively. In this paper we consider small deviations from the 
chemical equilibrium, Sfi <C T. In this case the quantities Ar m (m = 1,. . .,4) can be expanded in powers of 5/i and 
presented in the linear approximation as (see, e.g., (27l. l29l. l30l|): 

Ar m = X m Sfx, (12) 

where X m are the rates of the reactions (pQ)— <J3]) , some functions of the number densities and temperature. 



III. THE RELATIVISTIC HYDRODYNAMICS OF A SUPERFLUID NUCLEON-HYPERON MIXTURE 

In this and subsequent sections, the subscripts i and k refer to baryons (i,k = n,p, A, S). The summation is 
assumed over repeated baryon indices i and k. The subscript / refers to leptons (I = e, //); the subscript j runs over 
all particle species (J = n, p, A, E, e, fx); \i and v are the space-time indices. 

The relativistic hydrodynamics, describing a nucleon-hyperon mixture, composed of superfluid neutrons, protons, A 
and S _ hyperons, as well as normal electrons and muons, has been formulated in Ref. [271 ]. It is a direct generalization 
of the hydrodynamics of superfluid npe matter [16j . In this section we briefly discuss the main equations of this 
hydrodynamics. 
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A. The nondissipative hydrodynamics 



Let us discuss first the nondissipative hydrodynamics of supcrfluid nucleon-hyperon matter. Equations of superfiuid 
hydrodynamics include the continuity equations for particle species j 



where the particle four currents equal 



energy-momentum conservation law 



0. 



where the energy-momentum tensor T^" equals 

T i»> = (p + e ) u f u " + Pif v 



i k 



the second law of thermodynamics 

de = T dS + fa dn, + /i e dn e + /i M dn p 



(13) 

(14) 
(15) 

(16) 
(17) 



and a number of conditions for superfiuid components, which are specified below (see Eqs. (|20[) and (|2ip V In Eqs. 
(fl"3|) - (fT7|) 7]^ = diag(— 1, +1, +1, +1) is the special relativistic metric; rij is the number density of particle species j; 
e, T, and S are the energy density, temperature, and entropy density, respectively; P is the pressure, which is defined 
by the same formula as for nonsuperfluid matter [l5l flt| 



P = —e + fJ-irii + p e n e + n^rifj, + TS; 
u^ L is the four velocity of the normal (nonsuperfluid) liquid component, normalized so that 

= -1. 



(18) 



(19) 



(we assume that all the nonsuperfluid components move with the same velocity u^). Next, Yik is the symmetric 
relativistic entrainment matrix of nucleon-hyperon mixture, calculated in Ref. [321 ] for zero temperature and in Ref. 
[33l | for arbitrary temperature. The important property is that if some particles (e.g., neutrons) are nonsuperfluid 
then the related elements of this matrix vanish, Y n k = Ykn = (fc = n, p, A, £). In the nonrelativistic limit this 4x4 
matrix is expressed through the nonrelativistic entrainment matrix as [la . [l6| : — Pit / {iriiiTLk) , where m, is 
the mass of a free baryon species i (the matrix pa- is a generalization of the superfiuid density to the case of mixtures, 



see, e.g., Refs. [3J, |35j, |36j). The motion of superfiuid component of a species i is described by the four vector w 



which meets the condition 



(20) 



The potentiality of superfiuid motion is expressed as 



d" 



y£ } + q l A t * + p^ 



d» w v {£) + q t A v + (M u v 



(21) 



where qi is the electric charge of particle species i; A^ 1 is the four potential of the electromagnetic field. The potentiality 
condition (f2"Tj) is equivalent to a statement that there is a scalar function fa, satisfying (see Ref. fl6L l27l|) 



d»fa 



•jfa + qiA" + p i u fJ 



(22) 



The scalar fa is related to the wave function phase of the Cooper-pair condensate §i by the equality Vfa = V<I>i/2. 
In the nonrelativistic limit the spatial parts u and of the four vectors and transform into 



u 



W(i) = mi(V s (i) - V q ), 



(23) 
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where V q and V a U) = (V< fe — QjA)/ m,; are, respectively, the normal and superfluid velocities of the nonrelativistic 
theory of superfluid liquids [33, [38| . 

The hydrodynamics described above would be incomplete without an indication in what reference frame we define 
(measure) the main thermodynamic quantities (i.e., what frame is comoving). As was demonstrated in Ref. [l6j . the 
condition (|20j) dictates that the comoving is the frame where the four velocity equals = (1,0,0,0). In this 
frame, the basic thermodynamic quantities e, rij, and w^) (or Vifii) are defined by (see Eqs. (fT4|) . (fT6|) . and ([20]) ) 

fj = nj, (24) 
j t = Y ik w {k ) =Y ik (Vcf> k - q k A) , (25) 
T 00 = e. (26) 

All other thermodynamic quantities in nonequilibrium matter are the same functions of e, rij, and (or, equivalently, 
e, rij, and w^w^ k ^ L ) as in the full thermodynamic equilibrium. 

B. Viscous dissipation 

The main dissipative mechanisms in the pulsating nucleon-hyperon matter are the shear viscosity and effective bulk 
viscosity due to the nonequilibrium processes (H])-([4]). These are the two mechanisms which will be analyzed in this 
paper. 

The shear viscosity leads to an additional term in the energy-momentum tensor T^ v '. It now takes the form 



T" v = (P| £ )«V + Pif v 



+ Y ik 

where 



w (i) w (k) + Mi wf k) u v + fi k w v {l 



'sh ! 



(27) 



t^ = - v H^H vS (d 5 u 7 + d 7 u s -^ri lS d e uA; (28) 

jjhv = jj/ii/ _|_ u n u v j s ^j^g projection matrix; ?y is the shear viscosity coefficient. 

As has been already mentioned, the effect of the nonequilibrium processes dU)— (J4j) can be described in terms of 
the bulk viscosity formalism (for more details, see [13, H^, [3(| ) • However, for the analysis of sound modes it is more 
convenient to take these processes into account explicitly, by introducing sources into the right-hand sides of the 
continuity equations (|13[) . 

drffa = AS,, (29) 

similar to how it was done for normal (nonsuperfluid) npe matter in Ref. (39l | . Here ASj is a number of particles of a 
species j, generated in the unit volume per unit time in the reactions (|T|)— (J3J). Since we neglect the leptonic reactions, 
AS e = AS^ = 0. 

As a result, the hydrodynamics describing a viscous superfluid nucleon-hyperon mixture differs from the nondissi- 
pative hydrodynamics of Sec. IIIA only by the expression for the energy-momentum tensor (Eq. (f27|) instead of ([16])) 
and by the continuity equations (Eq. (|29p instead of (|13p ). Using this hydrodynamics, one can derive the entropy 
generation equation (see also p^.[27|) 

T dp (Su^) = 5fiAT- t% d^, (30) 

where we define 

Ar = Ari + Ar 2 + Ar 3 + Ar 4 = a Sfi, (31) 

with A = Ai + A 2 + A 3 + A 4 (see Eq. (fT2|) ). 

IV. SOUND WAVES IN THE ABSENCE OF DISSIPATION 
A. The main assumptions 

We assume that all the hydrodynamic velocities equal zero in the equilibrium, w' 1 = (1, 0, 0, 0) and = (0, 0, 0, 0) 
(i = n, p, A, E). (In principle, it is not necessary to have vanished all the components of the four vector It 
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is well known, that even in thermodynamic equilibrium a motion is possible with nonzero superfluid velocities [37| . 
This means that generally, the spatial components of the four vector can be nonzero. As for the time component 

w^j, it vanishes in any case, as it follows from Eq. ([20")) .) 

We assume also that in the perturbed matter the quasineutrality condition holds 

n p = n e + n,j, + ns, (32) 

in addition to the equilibrium condition © with respect to the fast reaction ([8]). Using Eq. (|32|) and the continuity 
equations H29p for protons, £~-hyperons, electrons and muons, as well as the fact that AS p — AS*,, — AS^ — AS^ = 0, 
one obtains [27| 

d^ Y pk w^ k) = df, Yz k W( k) ■ (33) 

In this paper we consider small deviations from the equilibrium state. Thus, we restrict ourselves to linear terms in 
perturbation. As follows from the normalization condition (| 19[) and Eq. (|20p . in the linear approximation the time 
components u° and w?^ in the perturbed matter remain the same, 

u° = 1, w° = 0, (34) 



while their derivatives vanish 



d^u° = 0, 0„ti$) = 0. (35) 
Using Eqs. (J34J) and ([35]), one gets from Eq. J33]) 



div [(Y pk - Yx k )w {k) ] = 0. (36) 

As we have already emphasized in Sec. IIIA, in the noncquilibrium matter any thermodynamic quantity (e.g., P 
or fjbj) is a function of the number densities rij, temperature T, and scalars w^W(k)^ (we choose T instead of e as an 
independent variable). In the linear approximation we can neglect the dependence on the quadratically small quantity 
w (i) w (k)v Moreover, in the strongly degenerate nucleon-hyperon matter, the dependence of P and \ij on T can also 



be neglected (see, e.g., [la,l22l,|4fj(). Consequently, P and fij are functions of only rij (j = n, p, A, S, e, fi). These six 
number densities are related by the conditions © and (|32[) . so that only 6 — 2 = 4 of them are independent. Thus, 
the pressure P and the chemical potentials jij depend on some four number densities (or their functions). As the 
independent variables it is convenient to choose the number density of baryons n& = n n + n p + n\ + ns, the number 
density of hyperons n# = n\ + n^, and the quantities n-£ n = n n + n n an d y — n e /n p (see Ref. ^3])- Notice that, 
in the thermodynamic equilibrium rib, nn, and y are additionally constrained by the two conditions (O and by 
the condition (jlip . In this case P and iij are functions of only one number density (e.g., rib). 



B. Equations governing the sound waves 



In this section we derive the system of equations describing sound waves in the superfluid nucleon-hyperon matter 
neglecting dissipation due to the nonequilibrium reactions ([TJ-© and shear viscosity. 

As we demonstrate below, the nonequilibrium reactions do not lead to dissipation in two limiting cases (i) either in 
the limit of slow reactions, when the total rate A of the reactions |l])-([4|) is negligible, so that they cannot change the 
matter composition during the pulsations excited by a sound wave (i.e. Ar m = 0, m = 1,. . .,4); (ii) or in the limit of 
fast reactions, when the reaction rate is so high that the pulsating matter is always in equilibrium with respect to the 
reactions so that the condition (jTTJ) is always satisfied. These are two cases which will be analyzed in what 

follows. 

The system of linearized hydrodynamic equations, describing sound waves, consists of (i) the condition (|36|) : (ii) 
momentum conservation law 

d t [(P Q +e Q )u + Li t0 Y tk w {k) ] = -VP, (37) 
following from Eq. (|15[k and (Hi) the four potentiality conditions for superfluid motion (for each baryon species) 

d t [/x n0 « + W(n)] = -V/x„, (38) 

d t [haoU + w w ] = -V^a, (39) 

d t [n P o u + w {p ) + q p A] = -V [fip + q p A°] , (40) 

d t [nvo u + I0 (s) + q^A] = -V + qzA°] . (41) 
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These conditions can be obtained from Eq. (j21j) with the help of Eq. ([35]) . The vector potential A and the scalar 
potential A can be excluded from Eqs. (|4T))) and (|41[) if one notices that = —q p and takes a sum of these formulas. 
Subtracting then Eqs. (|38|) and ([39| from the obtained sum and using the equilibrium condition (|9|), one gets 

#i [l»(2) + - 10(A) ~ W(n)] = 0. (42) 

In view of the condition ([6]), and the definition (fTO)) . it follows from Eqs. (|38|) and (f39|) that 

<9 t -t«(A)] = -VSfj,. (43) 

Eqs. ([HI) , and represent the five equations for five variables, w and W(j) (i = n, p, A, E). They 

completely describe the sound waves in the superfluid nucleon-hyperon matter under the condition that the pressure 
P{n b , njj, fisn, y)i neutron chemical potential n n {n b , njj, y), and chemical potential disbalance 5[i(n b , njj, "Em y) 
are known as functions of u and Wn\. Expanding these quantities in Taylor series near the equilibrium point and 
presenting P and u n as P = Pq + 8P and u n = /i„o + Sa n , one obtains 

ill' 

SP = — Sn b + 5n H + Sn Sn + — Sy, (44) 



dn 



b 



dn 



b 



dP 


■ Sn H - 


dP 


■ 5nz n - 


dP 


dn H 


dn-zn 


dy 


du n 


Sn H -+ 


du n 


Sns n -4 


dn a 


dn H 






dy 


dSfj, 
dn H 


Sn H - 


dSu, 


Sn^n - 


dSu, 
ay 



S^n = tj ~~ Sn b + 8n H + Sn^n + fy> (45) 



Su = - — Sn b + Sn H + ""^ Sn^ n + Sy, (46) 
dn b dn H dn Sn dy 

where in the last equation we take into account that S/j, — in the equilibrium state (see Eq. (fTTjl ). In Eqs. (144j) - (j4"6"|) 
Sn b , Suh, Sn-sn, and Sy are the deviations of the quantities n b , nu, iKn, and y from their equilibrium values n b Q, 
n Ho, n SnO? an d yo, respectively. In the Appendix these deviations are expressed through the velocities u and wu\ in 
the limit of slow and fast reactions. 

Assuming now, that the perturbations are harmonic (~ e M. u > t +' er ) } where ui is the pulsation frequency and k is the 
wave vector), the system of Eqs. (|36)) - (|38)) . (|42| . and (I43[) can be rewritten as 

(Y pk -Yj: k )w (k) =0, (47) 

iw [(P + e )u + n i0 Y lk w {k) ] = ik SP, (48) 

iw [u n0 u + W( n )] = -ikSu n , (49) 

W(A) + «>(»») = f(E) + t»(p), (50) 

[«)(„) -W(A)] = -ikSu. (51) 

In the limit of /asi reactions, when the condition fTTj) holds, Eq. ([ST]) can be further simplified, 

«>(n) - »0(A) = 0- (52) 



Eqs. (|47p - (f5Tj) together with (|44|) — (|46|) . as well as the expressions (|74j) ~(f76 |) for the case of slow reactions and 
(|77|) - (j79"l) for the case of fast reactions, allow one to determine the speeds of sound modes. 



C. Results for sound speeds 

From the analysis of Eqs. (J47|) -([5"l" j) it is clear, that the vectors u and Wu\ must be collinear with the wave vector 
A;: || k. Thus, the system of equations ([T7|) - ([5T]) has the form 

A-x = 0. (53) 

Here x is a vector that equals x = (u, u>(„), n>(p), U»(A), VJ(£)), where u = uk/k and = w^k/k; A is a 5 x 5 matrix, 
which elements depend on the thermodynamic quantities (and their derivatives) , on the relativistic entrainment matrix 
life, and on the frequency w and the wave number k. 

The system of equations (l53j) has a nontrivial solution if det A — 0. This condition results in a cubic equation on the 
speed of sound s = uj/k squared. Each of three solutions to this equation describes two sound waves, propagating with 
the same speed along and opposite to the wave vector k. In principle, these solutions can be written out analytically, 
but here we do not present them because they are too lengthy. 




FIG. 1: Speed of sound modes s in units of c versus T at m = Snto = 0.48 fm~ 3 for the third equation of state of Glendenning 
[4ll ]. Pulsation frequency is oj — 10 4 s _1 . Left panel shows three sound modes ('normal', 'sfl I', and 'sfl IF) in the limit of slow 
reactions. Right panel shows two sound modes ('normal' and 'sfl') in the limit of fast reactions. Baryon critical temperatures 
are indicated by vertical dots. Range of T where the limit of slow (left panel) and fast (right panel) reactions is invalid (for 
uj — 10 4 s _1 ), is shaded. 



In the case of slow reactions (when the condition \(d8u/ dun — dSfi/dnsn) *C u) holds, see Appendix), the three 
nonzero solutions to the cubic equation exist, corresponding to three different sound modes. In the case of fast 
reactions (when X(dS/i/dnH ~ dSfj,/dns n ) 3> uj), one of the three roots of the cubic equation vanishes, hence the 
actual number of sound modes is two (this conclusion can be also drawn from the fact that in the limit of fast 
reactions one can use Eq. ([52")) instead of (fSTj) , which does not depend on k and uj). 

The number of independent sound modes in these two limiting cases can be easily understood from the following 
reasoning. When all baryon species are superfluid, there are 5 velocity fields in the matter, u and Wuy In the limit 
of slow reactions the velocities are related by Eqs. (|4"T|) and (fB"0")) (these equations follow from the quasineutrality 
condition (|3"2")l and the condition of the equilibrium with respect to the fast reaction respectively). Thus, in this 
limit we have three independent velocities (5 — 2 = 3). Correspondingly, the number of independent sound modes is 
also three. In the limit of fast reactions, there is a constraint (|52[) in addition to the conditions (|4"T|) and ([50]) . This 
constraint is a consequence of the equilibrium condition (jlip with respect to the reactions Thus, the number 

of independent velocities (and sound modes) equals two (5 — 3 = 2). 

Figure [T] illustrates the results of numerical calculation of the sound speeds s (in units of c) as functions of temper- 
ature T in the limit of slow reactions (left panel) and fast reactions (right panel) . Here and below in this paper we use 
the third equation of state of Glendenning [4l[. The figure is plotted for the baryon number density n\> = 3nrj = 0.48 
fm -3 , where no = 0.16 fm -3 is the nucleon number density in atomic nuclei. The critical temperatures T C i (i = n, 
p, A, S) for transition of baryons to the superfluid state are poorly known. We take them to be T cn = 10 9 K, 
T cp = 2 x 10 9 K, T c a = 3 x 10 9 K, and T C £ = 6 x 10 9 K in accordance with some theoretical predictions (see, e.g., 
0) B EH) E3 ) • The relativistic entrainment matrix is taken from Ref. [33| (see also Ref. [Uj], where this matrix is 
calculated for T = 0). 

The shaded region in the left panel of Fig. [1] corresponds to temperatures, where the limit of slow reactions is 
invalid. Rather conventionally, we define it by the inequality X(dSfi/driH — dSfi/dn-£ n ) > uj/3. Similarly, the shaded 
region in the right panel of Fig. [T] indicates the range of temperatures, where the limit of fast reactions is invalid. We 
define this region by the inequality X(dSp,/driH — dSfj,/dnsn) < 3a;. 

Here and below we choose the frequency uj equal to u> = 10 4 s _1 . In fact, this value of uj is more appropriate for 
studies of pulsating neutron stars rather than sound waves, because it results in wavelengths of the order of stellar 
radius. However, we believe that the relatively simple analysis of sound waves with our choice of uj may provide some 
insight into the complex properties of global pulsations of neutron stars with superfluid nucleon-hyperon cores. 

The rates of the reactions Ai,. . .,A4, which constitute the quantity A = Ai + A2 + A3 + A4, are poorly known and 
depend essentially on the model of baryon-baryon interactions and on the many-body theory employed (see, e.g., 
[HI, HiJ [3^, [31], |42| ) . It is natural to expect that the rates A^ are order-of-magnitude comparable j28j. In our numerical 
calculations we, following Ref. [30| . take into account only the contribution to A from the reactions 0]) and ((2|) (i.e. 
we assume that A3 = A4 = 0). The rates A10 and A20 for nonsuperfluid matter are taken from Ref. [301 ] . The baryon 
superfluidity suppresses the reaction rates, which can be presented as Ai = Aioi?i and A2 = \2qR2- The reduction 
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factors i?i < 1 and R% < 1 are calculated using the formula (28) of Ref. [2{|. Though the authors of Ref. [29j 
proposed the formula for the reaction |T]), it remains valid for the reaction |(3J). In both cases the index i in that 
formula enumerate the reacting particles (i — n, n, p, E for the reaction ([1]) and i = n, p, p, A for the reaction ([2|)). 

As one can see from Fig. 1, at T < T cn (when all baryon species are superfluid) there are three sound modes in 
the limit of slow reactions and two in the limit of fast reactions. With further increasing temperature the number 
of sound modes decreases. As a result, at T > T c \ there is only one mode in both limiting cases. At T > T c s this 
mode becomes the ordinary sound in nonsuperfluid nucleon-hyperon matter (in fact, it transforms into the ordinary 
sound already at T > T c a, see the next paragraph). We term this sound mode normal; in the figure it is denoted as 
'normal'. Accordingly, the other modes are termed superfluid and denoted as 'sfl I', 'sfl II' (see Fig. 1, left panel) and 
'sfl' (see Fig. 1, right panel). Notice that, the speed of normal mode only weakly depends on T both in the limits 
of slow and fast reactions and practically coincides with the speed of sound of normal matter. On the contrary, the 
speeds of superfluid modes strongly depend on temperature and approach their asymptotic values only at T < 10 s K. 

Let us explain, how the number of sound modes changes with T in the limit of slow reactions (Fig. 1, left panel). 
Limit of fast reactions can be considered in a similar way. As we have already discussed above, at T < T cn there are 
five velocities, u and Wu\ (i — n, p, A, S). They are related by two conditions, (|47|) and (j50|) . Thus, the number of 
independent velocities (and sound modes) equals three. At T cn < T < T cp neutrons are normal, that is there are only 
four velocities in the system, u and wu> (i = p, A, £). These velocities are constrained by the only one condition (|47|) 
(the condition (|50[) is not applicable; it is valid only if all baryon species are superfluid). Thus, the number of sound 
modes remains equal three. Then, at T cp < T < T c \ and at T c \ < T < T c ^ a motion with, respectively, three (it, w\, 
and We) and two (u and w%) velocities is possible. These velocities are related by the condition (|T7|) . Consequently, 
there are two sound modes in the range T cp < T < T c \ and one in the range T c \ < T < T c j> In the latter case 
it follows from the condition (jl?)) that roj = 0. Hence, the hydrodynamic equations are formally the same as those 
for normal liquid, and the only sound mode coincides with the ordinary sound in nonsuperfluid matter. Finally, at 
T > T c £ all baryon species are normal and move with the same velocity u. Since this velocity is not constrained, 
there is (as it should be) only one sound mode. 



V. DAMPING OF SOUND WAVES 
A. Damping times: general equations 

In the previous section we analysed the sound modes in the limit of slow and fast reactions neglecting dissipation. 
In this section we allow for a weak dissipation in these limiting cases, which is due to the shear viscosity and 
nonequilibrium processes Our aim is the calculation of the characteristic damping times r of sound waves 

(the so called e- folding times). 

A few ways exist to calculate r. For instance, one can explicitly take into account the shear viscosity in Eqs. (147[) - (j5"Tj) 
and the next (complex) terms in the expansion of Snjj and Sn^ n into series in powers of X(d6fi/dnjj — dS^i/dn^ n )/uj in 
the limit of slow reactions and in powers of Lo/[\(d6fi/ dun — 06 fi/ dn-£ n )] in the limit of fast reactions (see Appendix). 
Then, solving the system of equations (|44)) - (|5T|) . one can find the small complex correction Ss to the speed of sound 
s, which is related to the damping time r by 



Another way to calculate r is to use the effective bulk viscosity formalism [27| , similar to how it was done in Ref. [16| 
for npe matter. 

We calculated r by both methods described above. However, here we present the most simple, third method of 
calculation. Of course, all three methods give the same results. 
Let us define the characteristic damping time as 



2E, 



puis 



(Epuh 



(55) 



where E pu \ s is the mechanical energy (per unit volume) of pulsations; (i^puis) is the rate of change of E pu \ s , averaged 
over the period 2tt/lo. Here and below angle brackets denote averaging over the pulsation period. The factor 2 in Eq. 
(|55[) appears because r is the e-folding time for the hydrodynamic velocities rather than for energy One can check 
that this definition of r coincides with that given by Eq. ([51]) . 

The mechanical energy of a pulsating superfluid matter can be easily found if we notice that -Epuis entirely transforms 
into its kinetic energy when the matter in the course of pulsations passes through the equilibrium point. Generally, 
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the kinetic energy -Ekin of the superfluid matter equals (we remind that ii<c and wu-j/fXi <C c) 

Skin = ^ {(P + e)u 2 + Y ik [m u)(fe)tt + ^kW(t)U + W(i)W( k )]} ■ (56) 

In the nonrelativistic limit and for two-component mixture this formula agrees with the expression (7) of Ref. [34| . 
Since in a sound wave the vectors u and Wu) are collinear with A; and can be presented as u = u a cos(a;A + kr) and 
W(i) = «J(j)a cos(wi + kr), one obtains for -E pu i s (u a = u a k/k and Wuy a = Wuy a k/k) 

Sp U i s = - {(P + e ) u 2 a + Y lk [fj, i0 w {k)a u a + {i k0 w (l)a u a + w^ a w^ a ] } . (57) 

Here we used the fact that at the equilibrium point the quantities u and W/{\ are equal to their amplitudes u a and 
I0(i) o , respectively. 

To find (Euuis) let us notice that all this energy goes into heat. Thus, using Eqs. (f5TJ)) and (pTLj) . one can write (see 
also Refs. [39|, l43j ) 

(Spms) = —(Tdp = -(\(Sn) 2 ) + (t£ d^uu). (58) 

It follows from this equation together with Eq. (155p that the characteristic damping time Tbuik due to the nonequilib- 
rium reactions (P)-(I1]) equals 

Tbuik Wf¥~r { ' 

while the characteristic damping time r s h due to the shear viscosity is 

2-Epuis / en \ 
T"sh = - , r, r. 60) 

As a consequence of Eq. (|T3[) (see Appendix) , in the limit of slow reactions 

/wc n2\ ^ /" ^A* t T 9S^i \ 2 

(A(0/i) ) = I — J&a + T; J//a + ^ JSna , (61) 

while in the limit of fast reactions 



(A(^) 2 )= . - ~ ax ~ ,J ^J ba + ^J Ha + ^J, na ) ■ (62) 
2s 2 X(dS t i/dn H - dSfi/dn^n) 2 \dn b dn H <9n s „ / 

In Eqs. (HU) and (JB2J) J fc a = mUa + ^YikW^a, Jh a = n H u a + YY, k w {k ) a + Y Ak w^ k ) a , and J En a = n s „u a + Y Efe W( fe)a + 
^nfe^(fe)a; s is the speed of sound calculated in the previous section neglecting dissipation; the factor 1/2 is a result 
of the averaging over the pulsation period. 

The dissipation rate of the mechanical energy due to the shear viscosity is the same in both limits, 

2 lu 2 

( T sh d H u ») = -^ J 1^2 "a- ( 63 ) 

To obtain this formula we used Eq. (f2"B"|) for rjt w . We see that (t^ duU v ) is formally given by the same expression as 
in the case of nonsuperfluid matter. 

As follows from Eqs. I|57[). (|58| . and (f6"Tj) - (|63p . E vvl \ s and (-Epuis) are the functions of the amplitudes u a and Wu\ a . 
Using the system of linear nondissipative equations (|44p -([5T j) . the amplitudes Wuy a can be expressed through u a and 
presented in the form 

W(i) a = ati(s) u a , (64) 

where cti(s) are some coefficients depending on s; they differ for each sound mode. In view of Eq. (|64p . one gets from 
Eqs. ([57]) . ([58]) . and (f6"l"j) - ([63"l) that i? pu i s ~ u 2 , and (E pu \ s ) ~ u 2 a . Hence, r is independent of u a (see Eq. (|55]l). 
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log 10 T (K) log 10 T (K) 

FIG. 2: The characteristic damping times r s h due to the shear viscosity versus T in the limit of slow reactions (left panel) 
and fast reactions (right panel). Dashed curves (marked with 'nfh') in both panels are obtained using the nonsuperfluid 
hydrodynamics, see the text. Other notations are the same as in Fig. 1. 



B. Results for damping times 

We numerically calculate the coefficients ati(s) from Eqs. (|44|) - (f5Tj) and thus determine the characteristic damping 
times r of sound modes in the limit of slow and fast reactions. Figures 2, 3, 4, and 5 illustrate the results of our 
calculations. These figures are plotted assuming the same microphysics input (the baryon number density, the critical 
baryon temperatures etc) as in Fig. 1. 

Figure 2 shows the characteristic damping times r s h due to the shear viscosity as a function of temperature T in the 
limit of slow reactions (left panel) and fast reactions (right panel). To calculate it is necessary to know the shear 
viscosity coefficient rj of superfluid nucleon-hyperon matter. This coefficient has not been considered in the literature 
so far. For definiteness, we take for r/ the shear viscosity of electrons and muons r\ eil — rj e + from Ref. [44J. In 
this reference it is shown that (for npe/x matter) the contribution of r/ efl to the total shear viscosity r\ is dominant. 
Notice, however, that this result was obtained under assumption that only protons are possibly superfluid (neutrons 
were treated as normal) . When plotting Fig. 2 we used the coefficient n efJi calculated for nonsuperfluid matter from 
Eq. (37) of Ref. The effects of baryon superfluidity on r s h are illustrated in the next figure. 

By the solid curves we show r s h calculated for each sound mode by means of Eq. (|60|) . Dashes demonstrate the char- 
acteristic damping times r n fh- s h calculated using the simplified model, the hydrodynamics of normal (nonsuperfluid) 
liquid. In this case there is only one mode in both limits; the corresponding curves are marked with 'nfh', which is the 
abbreviation of 'normal fluid hydrodynamics'. It is straightforward to demonstrate that r n fh- s h = 3(Po +£o)s 2 /(2w 2 f7) 
(see, e.g., Ref. [43|). As it should be, T s h for normal mode (in the figure it is marked 'normal') coincides with r n fh- s h 
at T > T cA (see Sec. IVC). 

It follows from Fig. 2 that r s h, calculated in the frame of relativistic hydrodynamics of superfluid mixtures, can 
strongly (by several orders of magnitude) differ from T n fh_ s h. It is interesting that the maximum deviation of T s h from 
Tnfh-sh at low temperatures (T < 3 x 10 8 K) is observed for normal mode (Fig. 2, left panel), though it is analogous 
to the usual sound in nonsuperfluid matter. 

At temperature T ks 5.55 x 10 s K the characteristic damping time r s h for one of the superfluid modes ('sfl I') 
becomes infinite. This is because at such temperature and for this mode the hydrodynamic motions occur in such 
a way that the normal component is always at rest (u a = 0), while the superfluid components pulsate around it. 
Mathematically, this means that the coefficients ati(s) in Eq. are infinite. It follows then from Eqs. (|57|) . (I60p . 
and (|63p that dissipation due to the shear viscosity is absent. 

Figure 3 presents the dependence of the characteristic damping times on T in the limit of slow reactions. Three 
panels correspond to three modes (from left to right: 'normal', 'sfl I', and 'sfl IF). The damping times Tbuik, calculated 
for each mode using Eq. ([55]) . are shown by solid curves. For comparison, by long dashes we show the characteristic 
damping times T n fh-buik due to the noncquilibrium reactions |l])-(|4]), which are obtained using the hydrodynamics of 
nonsuperfluid liquid. When plotting the long-dashed curve superfluidity of baryons was taken into account only at 
calculating the total reaction rate A. Because there is only one sound mode in the nonsuperfluid hydrodynamics, this 
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FIG. 3: The characteristic damping times r versus T for 'normal' mode (left panel), 'sfl I' mode (middle panel), and 'sfl IF 
mode (right panel) in the limit of slow reactions. Solid curves demonstrate Tbuik calculated from Eq. (|59[) : long-dashed curve 
(the same in all panels) describes the characteristic damping times T n fh_buik> calculated using the normal fluid hydrodynamics. 
Dot-dashed and short-dashed curves show damping times T sn i and r s h2 due to the shear viscosity (see the text for more details). 
Other notations are the same as in Figs. 1 and 2. 
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FIG. 4: The same as in Fig. 3 but in the limit of fast reactions. The damping times T sn 3> Tbuik and are not shown. 



curve is the same in all three panels. 

One can see from the figure that for the normal mode the solid and long-dashed curves practically (on the logarithmic 
scale) coincide; they differ by a factor of 2 or less. On the contrary, Tbuik for superfiuid modes can differ from T n fh-buik 
by orders of magnitude. 

Notice that, at T w 5.41 x 10 s K the damping time for the superfiuid mode 'sfl IF becomes infinite. The point is at 
such temperature and for this mode the condition 5fi — is always preserved during the pulsations. The denominator 
in Eq. (|59|) is then vanished and Tb u ik tends to infinity. 

By dot-dashes and short dashes in Fig. 3 we show, respectively, the characteristic damping times T s hi and T sn 2 due 
to the shear viscosity. The times T s hi are the same as in Fig. 2 (left panel); they are plotted for the shear viscosity 

V = r hn °f nonsuperfluid matter (see equation (37) of Ref. [44j). To obtain the dependence T s h2(T') we assume that 

V = Vefj, as before, but additionally take into account the reduction of rj efi by proton superfluidity The reduction 
factor was calculated from Eq. (83) of Ref. [44j. It is worth noting that this formula is obtained for nucleon npe/j, 
matter and does not imply the superfluidity of other baryon species except for protons. Thus, the dependence T sn 2(r) 
only qualitatively describes possible effect of baryon superfluidity on T sn . 

It follows from the analysis of Fig. 3, that at high enough temperatures T > 3 x 10 8 K, the dissipation due to 
the shear viscosity is negligible in comparison to that due to the nonequilibrium processes dU) — - Moreover, the 
threshold density, at which Tbuik ~ T s h, only weakly depends on rj. 
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8.5 9 9.5 10 

log l0 T (K) 

FIG. 5: Characteristic damping times rbuik versus T in the limit of slow (T < 1.4 x 10 9 K) and fast (T > 2.3 x 10 9 K) reactions. 
The region of T where the rate of nonequilibrium reactions is intermediate (i.e., where both limits are invalid), is shaded. 

Figure 4 is analogous to Fig. 3 but is plotted for the limit of fast reactions. The solid curve in the left panel 
demonstrates the dependence Tbuik(T') for the normal mode ('normal') and in the right panel - for the supcrfiuid 
mode ('srl'). The long-dashed curves are the same in both panels and describe the damping times t^^^^T) . In the 
limit of fast reactions the damping times due to the shear viscosity are much greater than Tb u ik and are not shown. 

As in the limit of slow reactions, Tbuik for the normal mode does not differ substantially from T n f n _buik (no more 
than by a factor of 2). However, at high enough temperatures Tbuik for superfluid mode is two orders of magnitude 
smaller than T n fh-buik- 

The analysis of Fig. 4 shows that the damping time for superfluid mode becomes infinite at the same temperature 
T « 5.41 x 10 8 K as in the limit of slow reactions (the corresponding peak is in the shaded region where the limit of 
fast reactions is not applicable for the chosen pulsation frequency u = 10 s _1 ). It is easy to understand why this 
peak fall on the same temperature. As we already discussed above, the infinite damping time Tbuik means that the 
condition Sfi — holds in the pulsating matter (at given temperature). In this case, the nondissipative Eqs. (|47 p - (|5"Tj) 
in the limits of slow and fast reactions coincide. Thus, the functional dependence of the amplitudes W(i) a on u a is the 
same in both limits (see Sec. VA). Consequently, as follows from Eqs. (jrTTj) and (frj2"|) . if = in one limit, then it 
vanishes in the other limit. 

Fig. 5 is a combined plot, showing the characteristic damping times Tbuik as functions of T for the limits of slow 
and fast reactions. The shaded region (at T ~ 2 x 10 9 K) corresponds to intermediate temperatures at which both 
limits are inapplicable. In this region the pulsation energy dissipates on a time scale of order of the pulsation period. 
Notice that, at T > 3 x 10 8 K the characteristic damping time of pulsations r = Tbuik + T s h practically coincides with 
Tbuik, since at such temperatures the contribution of the shear viscosity to dissipation is negligible. 

VI. SUMMARY 

In Ref. [27J the relativistic dissipative hydrodynamics was suggested to describe superfluid nucleon-hyperon matter 
in the cores of massive neutron stars. Using this hydrodynamics, we analyse the sound waves, which are the simplest 
example of pulsations in such matter. 

We demonstrate that in the limit of slow reactions (H])-(|l]) (when the composition of pulsating matter is practically 
unaffected by these reactions) there are three sound modes: one normal and two superfluid. In the opposite limit 
of fast reactions (when the pulsating matter is nearly at equilibrium with respect to the reactions ([T])-6])), only two 
sound modes exist: the normal one and the superfluid one. In the intermediate case the sound waves cannot propagate 
because they are damped on a time scale of order of the pulsation period. 

The speed of normal sound mode in both limits is practically independent of temperature and coincides with the 
sound speed for nonsuperfluid matter. This mode turns into the ordinary sound at T > T c \. On the contrary, the 
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speeds of superfiuid modes strongly depend on temperature and vanish before the transition of matter to the normal 
state. 

We analyse also the characteristic damping times r of sound modes (Figs. 2, 3, 4, and 5). We allow for the two 
main dissipative mechanisms: damping due to the shear viscosity and due to the nonequilibrium reactions (HJ) — d3I) 
(these reactions generate the effective bulk viscosity). We demonstrate that (i) the damping times t for normal and 
superfiuid modes can differ from each other by orders of magnitude; (ii) at T > 3 x 10 8 K the damping due to 
the nonequilibrium reactions (HJ) — (J3J) is the dominant mechanism of dissipation; this result is nearly insensitive to an 
actual value of the shear viscosity coefficient rj. 

In addition, we compare r with the damping time T n fh, calculated using the ordinary nonsuperfluid hydrodynamics, 
but taking into account the effects of superfluidity on the shear viscosity and on the rates of the reactions ©-([I]). 
We show that (Hi) r approximately (up to a factor of two) coincides with r n fh only for the normal mode and under 
the condition that the shear viscosity can be neglected (i.e. T > 3 x 10 s K). In other cases (for the superfiuid modes 
and for the normal mode at T < 3 x 10 8 K) t can differ from r n fh by several orders of magnitude. 

The results listed above are obtained from the analysis of sound waves in the superfiuid nucleon-hyperon matter. 
However, they can serve as an indication that the effects, related to difference between the superfiuid and normal 
fluid hydrodynamics, can also be very important in studies of global pulsations of superfiuid neutron stars, essentially 
modifying their damping times. 



Appendix 

Let us calculate the quantities Snb = rib — ri&o, Srin — Uh — fiHO, Sn^ n = nsn — Tis n o, and Sy — y — yo, entering 
Eqs. (3U) - ([15} for SP, S[i n , and <5/x, respectively. For that, we make use of the continuity equations (|2T))) . assuming the 
perturbations are harmonic. In the linear approximation the continuity equations for leptons (electrons and muons) 
have the form (I — e, jS) 

iuSm + ikniou = 0, (65) 

where we put AS/ = 0, because the leptonic reactions are slow. From these equations it follows that 

Sy = 0. (66) 

The continuity equations for baryons, hyperons, and S _ -hyperons with neutrons can also be obtained from Eq. (|29[) 

iu> 5rib + ikJb — 0, (67) 
ito Stih + ikJ h = A5/i, (68) 
iuj Sns n + ikJ-£ n — —XS/i. (69) 

Here we used Eq. (|3"TT) . and introduced the notations Jb = ribU + J2i YikW(k)) Jh = nnu + Y^kW(k) + Y\kW(k)'i an d 
Jsn = ns„u + Y^kW(k) + YnkW(k)- Solving now the system (|67)) - (|69|) taking into account Eqs. (|46|) and ([66]) . one gets 
for 5rib, dn,H, and 5n^ n 

Sn b = (70) 

Sn ^ _ iujJh + A[J b d6n/drib + (J H + Jsn) dSfi/dn^ n ] 

lu \(d5fi/driH — 85 /j,/ drtYto) — ito 

_ k iujJxn - A[J fc d5ti/dn b + (Jg + J Sra ) dS[i/dn H } 

lj \(dS[i/dnH — dSn/dn-z n ) — ioj 

Using these equalities, Eq. (I46|) for 5/i can be rewritten as 

s = k / ^J b + ^-J H + d5 ^ J z ^ (73) 

In the limit of slow reactions, when the total rate A of the reactions CD)-([I]) is small, that is X(ddfi/dnji — 
dSfi/dn-^n) <C u>, one has 

Sn b = -Ub/uJ, (74) 

5tih = —kJii/uj, (75) 

Sn^n = -kJ^ n /uj. (76) 
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In the limit of fast reactions, when \(d6fi/driH — dSfi/dn^ n ) 3> w, one obtains 



6n b 



kJ b 



(77) 



k(J b dS^i/dut, + (J H + Js,i) dS^/dnxn) 
uj(dS{i/dn,H — dSfi/dnsn) 
k(J b dSfi/dn b + (Jjj + J S n) dSfi/dn H ) 
uj(dS[i/driH — dd/s/dnsn) 



(78) 



(79) 



One sees, that in the both limits Sn b , <5n#, and Sns„ are real- valued and do not depend on the total reaction rate A. 
Correspondingly, the sound speeds are also real- valued, or, in other words, the dissipation is absent. The dissipation 
due to the weak nonequilibrium processes ©-(HI) can be taken into account by considering the next (complex) terms 
in the expansion of Suh and 5n^ n into series in powers of X(ddn/dnj] — dSfi/dn^ n )/u; in the case of slow reactions 
and in powers of uj /[X(dSfj,/ 'dnn — dS^i/dns n )] in the case of fast reactions. 
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